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Abstract. We introduce a novel migration process, the target process. This process 
is dual to the zero-range process (ZRP) in the sense that, while for the ZRP the rate 
of transfer of a particle only depends on the occupation of the departure site, it only 
depends on the occupation of the arrival site for the target process. More precisely, 
duality associates to a given ZRP a unique target process, and vice-versa. If the 
dynamics is symmetric, i.e., in the absence of a bias, both processes have the same 
stationary-state product measure. In this work we focus our interest on the situation 
where the latter measure exhibits a continuous condensation transition at some finite 
critical density p c , irrespective of the dimensionality. The novelty comes from the 
case of asymmetric dynamics, where the target process has a nontrivial fluctuating 
stationary state, whose characteristics depend on the dimensionality. In one dimension, 
the system remains homogeneous at any finite density. An alternating scenario however 
prevails in the high-density regime: typical configurations consist of long alternating 
sequences of highly occupied and less occupied sites. The local density of the latter is 
equal to p c and their occupation distribution is critical. In dimension two and above, 
the asymmetric target process exhibits a phase transition at a threshold density po 
much larger than p c . The system is homogeneous at any density below p$, whereas 
for higher densities it exhibits an extended condensate elongated along the direction 
of the mean current, on top of a critical background with density p c . 
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1. Introduction 

In recent years many studies have been devoted to nonequilibrium statistical-mechanical 
models yielding condensation, such as zero-range processes (ZRP) [1, 2, 3, 4, 5, 6, 7], 
dynamical urn models [8, 9, 10, 11, 12], and mass transport models [13]. In all these 
models the condensate manifests itself by the macroscopic occupation of a single site by 
a finite fraction of the whole available mass. 

The ZRP is the simplest of these models. It is a driven diffusive system with 
multiple occupations, such that particles hop from site to site on a lattice, with a rate 
which only depends on the occupation of the departure site. The distribution of the 
particles among the sites in the stationary state is given by a product measure, which 
is explicitly known in terms of the rate defining the model, irrespective of the geometry 
of the system and of the asymmetry of the dynamics [6, 14, 15]. This property of 
the stationary-state measure of the ZRP favors the condensation phenomenon. The 
product structure indeed prevents the condensate from being shared by more than one 
site [5, 7, 13]. 

Dynamical urn models, also called migration processes in the probabilistic litera- 
ture [16], can be viewed as generalizations of the ZRP, where the rate at which a particle 
is transferred from a departure site to an arrival site now depends on the occupations 
of both sites. Migration processes do not have a stationary-state product measure in 
general. 

In the present work we introduce a special migration process, the target process. 
This process is novel to our knowledge. It has no stationary-state product measure, 
except in the case of a symmetric dynamics, i.e., in the absence of a bias, leading to an 
equilibrium stationary state. The class of target processes is dual to the class of ZRP, 
in the sense that the roles of the departure and arrival sites are interchanged: the rate 
basically depends on the occupation of the departure (source) site for a ZRP, and on 
the occupation of the arrival (target) site for a target process. More precisely, to a given 
ZRP is associated by duality a unique target process, and vice-versa. In the case of a 
symmetric dynamics, these two dual processes have the same stationary-state product 
measure. 

Our aim is to study the structure of the nonequilibrium stationary state of the 
asymmetric target process, and especially the fate of the condensation phenomenon. In 
Section 2 we give a detailed definition of the class of target processes. We then focus our 
attention onto the particular target process which is dual to the 'canonical ZRP' studied 
e.g. in [2, 3, 4, 5]. In the case of symmetric dynamics, the target process thus constructed 
has the same stationary-state product measure as the ZRP, and therefore the same 
continuous condensation transition. In the asymmetric case, however, the stationary- 
state measure of the target process is not a product measure, and exhibits non-trivial 
correlations in general. This basic difference between the ZRP and the target process 
manifests itself more drastically at high density and in low dimensionality. Sections 3 
and 4 are devoted to a thorough theoretical and numerical study of the one-dimensional 
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situation. The system remains homogeneous at any finite density, and presents an 
alternating structure which is more and more pronounced as the density is increased. 
The two-dimensional asymmetric target process on the square lattice is the subject of 
Section 5. It exhibits an unconventional type of condensation phenomenon, with a 
transition at a threshold density po much larger than p c , and an extended condensate 
elongated along the direction of the bias. An analogous scenario presumably generically 
holds on other lattices and in higher dimension as well. Section 6 contains a Discussion. 

2. The target process 

2.1. Migration processes and urn models: a reminder 

The definition of a migration process appeared first in the probabilistic literature (see 
e.g. [16]). For the time being, we restrict the discussion to the one-dimensional case. 
Consider a system of iV particles distributed among M sites, with periodic boundary 
conditions. Let N m be the number of particles at site m = 1, . . . , M. A migration 
process (or dynamical urn model) is defined by the following Markovian dynamics: (i) 
a departure (source) site d is chosen at random, uniformly among the M sites of the 
system; (ii) an arrival (target) site a is chosen among the neighbors of d. To be specific, 
the right neighbor (a — d + 1) is chosen with probability p, whereas the left neighbor 
(a = d — 1) is chosen with the complementary probability q — 1 — p; (iii) a particle is 
transferred from site d to site a at a rate W k ,i which only depends on the occupations 
k = N d and I = N a of the two sites involved. Of course, one has 

Wo,i = 0, (2.1) 

since no particle can be removed from an empty site. 

A one-dimensional migration process is therefore entirely defined by the bias p and 
the rate Wk,i- A natural question is the following: What are the conditions on the bias p 
and the rate W^i for a one- dimensional migration process to have a stationary-state 
product measure? This question was first addressed in [17] (for a review, see [12]). 

In this context, a stationary-state product measure means that the probability of 
any configuration of the system in its stationary state has the form 

P(N U ...,N M ) = — L_ p Nl ... PNm S(N 1 + --- + N m , TV), (2.2) 

where the factors pk are arbitrary, and the partition function Zm,n reads 

Zm,n= E Pn 1 ..-Pn m S(N 1 + --- + N m ,N). (2.3) 

N U ...,N M 

The answer to the above question is as follows: 

• For symmetric dynamics, i.e., when p = 1/2, the stationary state has a product 
measure characterized by the factor pk, if the rate Wk,i obeys the condition 

Pk+iPiW k+h i = p k pi + iW l+hk . (2.4) 
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The resulting stationary state is an equilibrium state. The relation (2.4) expresses 
detailed balance with respect to this equilibrium state. 

• For asymmetric dynamics, i.e., when p ^ 1/2, the stationary-state has a product 
measure if the rate W k j obeys both (2.4) and the following condition: 

W Kl - W kfl = W t , k - W lfl . (2.5) 

The stationary state is a genuine nonequilibrium steady state. The condition (2.4) 
does not express detailed balance any longer, albeit a weaker property [12], referred 
to as pairwise balance [18]. 

The partition function Z M>N associated with the product measure (2.2) can be 
rewritten, using an integral representation of the Kronecker delta function, as 

^Wd^W. (2 - 6) 

where the generating series of the factors p k reads 

P(z) = J2PkZ h . (2.7) 

k>0 

The product measure (2.2) is therefore characterized by M, N, and the factor p k 
or, equivalently, by the generating function P(z). For a homogeneous system in the 
thermodynamic limit, where M and iV are simultaneously large, with a fixed density 
p = N/M of particles per site, the contour integral in (2.6) can be evaluated by the 
saddle-point method. The saddle-point value z, which is to be identified with the 
fugacity in the grand canonical ensemble, is related to the density p through the equation 

The distribution f k = Prob{iV\ = k} of the occupations of any given site of the system 
can be derived by summing the probability (2.2) over N 2 , . . . , N M . We thus obtain 

, Z M -i,N-k (r. n \ 

fk = Pk — = • (2.9) 

In the thermodynamic limit, this expression simplifies to 

A - fg- (2.10) 
2.2. The example of the zero-range process 

In the present context, the ZRP appears as the special case of a migration process where 
the rate W k ,i only depends on the occupation of the departure site: 

W k ,i = u k , (2.11) 

with u = 0, by virtue of (2.1). The condition (2.5) is then automatically satisfied, 
irrespective of the bias p. Equation (2.4) yields the following relation between the 
rate u k and the factor p k : 

Pk = Ufc+iPfc+i, (2.12) 
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up to a multiplicative constant, which we set equal to unity by an appropriate choice 
of time unit. The corresponding factor p k can be expressed in terms of the rate u k as 
follows: 

Po = l, Pk = (k>l). (2.13) 

Ui...U k 

Reciprocally, to a given stationary-state product measure characterized by the factor p k , 
there corresponds a unique ZRP dynamics (up to a choice of time unit), whose rate reads 

u k = 5*=1 (k > 1). (2.14) 

Pk 

2.3. Characterization of migration processes with stationary- state product measure 

We now give an explicit characterization of the migration processes which admit a 
stationary-state product measure. 

For symmetric dynamics, and for a given factor p k , the most general form of the 
rate obeying (2.4) reads 

W*^— S*-i,J (k>l), (2.15) 
Pk 

where the ratio p k _\jp k is nothing but the rate u k of the corresponding ZRP, given 
by (2.14), whereas S k j is a symmetric function of k and 

S k ,i = S l>k . (2.16) 

Besides the factor p k , the rate W k j depends on an arbitrary symmetric function S k j of 
two indices. 

For asymmetric dynamics, and for a given factor p k , the most general solution 
of (2.4) and (2.5) is determined by the one-dimensional array of rates a k = W kt o [19]. 
Note that a k is the rate at which an empty site is refilled, by receiving one particle from 
a non-empty neighboring site containing k > 1 particles, and that (2.1) implies ao = 0. 
The above property can be shown as follows. Consider all the indices k and I for a fixed 
value of the sum k + I = n, and introduce the quantities 

A k = p k Pn-kW k , n - k (k = 0,...,n). (2.17) 
Equations (2.4) and (2.5) respectively become 

A k = A n+1 _ k , A k - p k p n - k a k = A n _ k - p k p n - k a n _ k . (2.18) 
Combining these two equations yields 

A k+ i - A k = p k p n -k(a n -k - a k ). (2.19) 
The solution of this inhomogeneous difference equation with initial value A = reads 

k 

A k = ^ Pk-mPn-k+mi^n-k+m ~ Ctk-rn) ■ (2.20) 
m=l 

We are thus left with the following expression for the rate W k y. 

1 k 

W k l = Pk-mPl+m(0li +m - ttfc-m)- (2.21) 

PkPl m=l 

Besides the factor p k , the rate W kt i depends on an arbitrary function a k of one index. 
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2.4- Definition of the target process 

We define the target process as the migration process where the rate 

W k;l = (1 - S k , )v t (2.22) 

essentially depends on the occupation I of the arrival (target) site. Note the dual 
character of this definition with respect to the definition of the ZRP, in that the roles 
of the departure and arrival sites have been interchanged. There is, however, a key 
difference between the two models, coming from the presence of the constraint (2.1). 
The latter, which explicitly enters (2.22) through the factor (1 — 5 k:0 ), implies that the 
rate of the target process actually also bears some dependence on the occupation k 
of the departure site, as it is constrained to vanish if k — 0. For the ZRP the same 
constraint does not change the fact that the rate u k only depends on the departure site; 
it just imposes uq — 0. 

For symmetric dynamics, the target process always has a stationary-state product 
measure, for any choice of the rate v k . Equation (2.4) yields the following relation 
between the rate v k and the factor p k : 

Pk+i = v k p k , (2.23) 

up to a multiplicative constant, which is again set to unity by an appropriate choice of 
time unit. The factor p k of the stationary-state measure can thus be expressed in terms 
of the rate v k as follows: 

p = l, p k = v ...v k - 1 (k>l). (2.24) 

Reciprocally, to a given stationary-state product measure, characterized by the factor p k , 
there corresponds a unique symmetric target process (up to a choice of time unit), whose 
rate reads 

v k = ^±1. (2.25) 
Pk 

Equations (2.12) and (2.23) show that the ZRP with rate u k and the target process 
with rate v k have the same stationary-state product measure, i.e., the same factor p k , if 
the rates obey 

v k = — , (2.26) 
u k+ i 

up to a multiplicative constant. A target process and a ZRP related by this condition 
are hereafter named dual to each other. 

It is interesting to consider a more general class of migration processes, where the 
rate has the form 

W^ = (l-Mw (2-27) 

For symmetric dynamics, this model again has a stationary-state product measure for 
any choice of u k and v\. Equation (2.4) indeed yields 

u k+1 p k+1 = v k p k , (2.28) 
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up to a multiplicative constant. The factor p k can thus be expressed in terms of the u k 
and v k as follows: 

Pd = l, p k = V °--- Vk - 1 (k>l). (2.29) 
u x ...u k 

This expression shows that the factor p k only depends on the ratio Vk/uk+i- This 
class of processes interpolates between the ZRP and the target process, which are 
respectively recovered as the special cases where v\ — 1 and u k — 1. We finally notice 
that the rate (2.22) of the target process, and more generally the rate (2.27), is of the 
form (2.15), with 

S k ,i = v k v h (2.30) 

For asymmetric dynamics, the condition (2.5) is very stringent. The target process 
has a stationary-state product measure if and only if the rate v k only assumes two values, 
according to whether k is zero or not: 

Vk = \v (k>l). (2 - 31) 

This also holds for the more general process defined by the rate (2.27). 

For a generic asymmetric target process, where the rate v k is not of the form (2.31), 
the stationary-state measure is not a product measure. It is not known explicitly, and 
can be expected to be a non-trivial correlated measure in general. 

2. 5. The 'canonical target process ' considered in this work 

The condensation phenomenon in the ZRP is usually investigated using the rate [2, 3, 
4, 5]: 

= ! + £ (A;>1), (2.32) 

where the control parameter b is a measure of the strength of interactions. The 
minimality and exemplarity of this choice of rate suggest to call this model the 'canonical 
ZRP for condensation', or 'canonical ZRP', for short. 

Throughout the following we focus our attention onto the target process dual to 
the canonical ZRP. Its rate v k is therefore related to the rate (2.32) by the duality 
relation (2.26). We thus obtain 

*=*TTTT ( * £0) ' (2 ' 33) 

We name this process the 'canonical target process'. 

For symmetric dynamics, it has already been shown above that the target process 
has the same stationary-state product measure as the ZRP. It therefore exhibits the 
same condensation transition. We now give a brief reminder of the properties of this 
stationary-state measure. 
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In the absence of interactions (6 = 0), the rate reads uj, = Vk = 1, irrespective of k. 
We have therefore pk = 1, so that P(z) = 1/(1 — z). The fugacity z and the density p 
are related by 

z = ^+T' ' = T^7 < 2 ' 34 > 

The distribution of the occupations (2.10) is a geometric distribution: 

f k = (1 - (2.35) 

In the situation where 6 is positive, the rate u k is a decreasing function of the 
occupation k, so that particles hop less easily out of more occupied sites. Accordingly, 
the rate Vk is an increasing function of k, hence particles hop preferentially towards more 
occupied sites. The rate (2.32) or (2.33) therefore corresponds to attractive interactions 
between particles. The model exhibits a trend toward segregation, which leads to a 
thermodynamical condensation transition if b is strong enough. It turns out that many 
characteristics of this condensation transition are universal [1, 2, 3, 4, 5, 6, 7]: they only 
depend on the asymptotic behavior of the rate u k (or Vk) at large k, i.e., essentially on 
the value of b. 

For the choice of rate (2.32) or (2.33), the factor pk of the stationary-state product 
measure (see (2.13) or (2.24)) reads 

pk = J f f + P k l =b f\l- u) b ^u k du, (2.36) 

so that 



T(k + b+l) 

P(z) = b / ^ - 1 — du. (2.37) 

K ' k \ - zu v 1 

The factor pk falls off as a power law with exponent b at large k: 

Vk » (2.38) 

The canonical ZRP has a condensation transition in the thermodynamic limit whenever 
the first moment of the factor pk, 

P'(l) 

Pc = E*P*= pTTV. ( 2 - 39 ) 

k>i 

is convergent. We have P(l) = b/{b - 1) and P'(l) = b/((b - 1)(6 - 2)), so that the 
critical density p c is finite for b > 2, and reads 

* = S§ = rh- (2 - 40) 

This critical density separates a fluid phase and a condensed phase: 

• At the critical density (p = p c , z = 1), the occupation distribution reads 
x - Pk b - 1 

fk =Pjl)=— pk (2 ' 41) 
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(see (2.10)). In particular the fraction of empty sites is 

fo = (2.42) 
whereas the distribution falls off as a power law for large occupations: 

The statics and the dynamics of the model exhibit many features of critical 
phenomena, including scaling and universality. 

• In the fluid phase (p < p c , z < 1), the occupation distribution fk falls off 
exponentially. 

• In the condensed phase (p > p c ), for a large finite system, the particles are arranged 
so as to form a uniform critical background and a macroscopic condensate, typically 
occupying one single site and consisting of N — Mp c = M(p — p c ) excess particles. 

The stationary-state measure of the asymmetric canonical target process is not 
a product measure in the presence of interactions, i.e., for any non-zero value of the 
parameter b, because the rate (2.33) is not of the form (2.31). This absence of a product 
measure also holds for the asymmetric target process on higher-dimensional lattices. The 
derivation of the conditions (2.4), (2.5) given in [12] could indeed easily be extended to 
the case of a biased one-particle dynamics on any higher-dimensional lattice. The rest 
of this paper is devoted to a detailed investigation of this model. The study of the one- 
dimensional totally asymmetric model is presented in Sections 3 and 4, whereas Section 5 
is devoted to a maximally asymmetric form of the model on the two-dimensional square 
lattice. 



3. One-dimensional target process: theoretical analysis 

In this section, we consider the target process defined by the rate (2.33) in the totally 
asymmetric one-dimensional case (p = 1). The situation of most physical interest will 
turn out to be the regime of a high density. 

We begin the analysis of the model by exploring the consequences of the existence of 
a conserved current. Consider a large finite system, with periodic boundary conditions. 
In the stationary state, the mean current of particles through the system is conserved: it 
assumes the same value J through every bond. The current through the bond between 
sites m and m + 1 is given by the mean value of the corresponding rate: 

■J=(W Nm , Nm+1 ). (3.1) 

The existence of this conserved current is expected to ensure some degree of homogeneity 
of the stationary state. 

To illustrate the method, let us first consider the totally asymmetric ZRP, with 
rate (2.32). The current therefore reads 

■J=(unJ. (3.2) 
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Consider a typical configuration in the condensed phase (p > p c ). The site where the 
condensate is located contains a macroscopic number of particles, so that (3.2) yields 
J = 1 for the bond to the right of the condensate, up to a negligible finite-size correction 
of order 1/M. For all the other bonds, the relations (2.10) and (2.12) yield 

J =J2 u kfk =z, (3.3) 
k>i ^~y~^ 

Zjk-l 

where z is the fugacity. Equating the two above expressions for the current, we recover 
well-known results for the totally asymmetric ZRP, i.e., z = z c = 1 and J = 1 throughout 
the condensed phase. 

Let us now turn to the totally asymmetric target process, with rate (2.33). 
Equation (3.1) for the current now reads 

J=((l-SN m ,o)v Nm+1 ). (3.4) 

At variance with (3.2), this expression involves the occupations of two consecutive sites. 

The high-density regime turns out to be the situation of most physical interest. In 
this regime, at least some of the sites must have large occupations. We are therefore led 
to distinguish between two types of sites: 

• B-sites (B for big), whose occupation is large, of the order of p. 

• S-sites (S for small), whose occupation is small and fluctuating. 

This distinction will be kept at a heuristic level throughout the following. A typical 
high-density configuration therefore consists of four types of bonds: BB, BS, SB, and SS. 
The typical values of the current in each type of bond obey the following inequalities in 
the high-density limit: 

Jss < f B < Jbb- (3.5) 

^BS 

Each factor in (3.4) is indeed less than unity, and approaches unity in the limit where 
the involved occupation goes to infinity. In the high-density limit we have therefore 
Jbb = 1, whereas Jbs = (vN m ) and Jsb = (1 — <^v m ,o), where m is a typical S-site, 
and Jss is smaller than the last two expressions. 

Thanks to the inequalities (3.5) the presence of a stable isolated condensate is 
excluded. This would indeed correspond to a current profile where the two bonds on 
either side of the condensate carry currents J = Jsb and J = Jbs which are significantly 
higher than the background current J = Jss of all the other bonds. This non-uniform 
current distribution with a point defect would have the effect that the condensate 
would soon dissolve into the background. The inequalities (3.5) actually only leave 
out two scenarios for typical stationary-state configurations in the high-density limit. 
Both scenarios, described below and illustrated in Figure 1, correspond to spatially 
homogeneous phases. We are therefore led to predict that the asymmetric target process 
has no condensation transition at any finite density in one dimension. 
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Figure 1. Typical configurations in the two possible scenarios for the stationary 
state of the asymmetric one-dimensional target process at high density. Left: 
uniform configuration. Right: alternating configuration. 



Uniform scenario. In the uniform scenario, shown in the left panel of Figure 1, 
typical configurations are entirely made of B-sites, whose typical occupation is 
around the density p. These configurations carry a current J = Jbb, he., 

J=l (3.6) 

in the high-density limit, up to a correction of order 1/p. 

Alternating scenario. In the alternating scenario, shown in the right panel 
of Figure 1, typical configurations are alternating, i.e., they have the form 
BSBSBSBS... 

The asymptotic value of the current J = Jbs = <^sb through a perfect alternating 
structure in the high-density limit can be evaluated as follows. Most particles belong 
to B-sites, whose typical occupation is very high, around 2p. Every B-site therefore 
acts as a reservoir, so that the occupations of the S-sites evolve independently from 
each other. The master equations for the occupation distribution fk{t) of any given 
S-site read 

^jj = fk+i + Ufc-i/fc-i - (1 + v k )fk (k > 1), 

^ = fi-v fo. (3.7) 

The stationary-state solution of these equations is such that fk+i = Vkfk- It is 
therefore proportional to the factor pk- The properly normalized solution is given 
by (2.41). To sum up, the occupations of the S-sites are independent variables, 
whose distribution coincides with the critical occupation distribution of the dual 
ZRP or, equivalently, of the symmetric target process. We have 

J=Y. v kfk = l-k- (3.8) 

Jk+l 

Using (2.42), the current in the high-density limit reads 

J = \. (3.9) 
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This alternating scenario holds a priori whenever b > 1. The latter condition 
corresponds to P(l) being finite, so that (2.41) is a properly normalized probability 
distribution. 

In order to have a hint on which of the two above scenarios is preferred, it is 
interesting to first consider the simple example of a system of two sites. In this case, 
there is only one kind of move for the particles, namely from one site to the other, 
and only one degree of freedom, the occupation Ni = k of site number 1. Indeed 
the occupation of site number 2 reads N 2 = N — k. The stationary-state occupation 
distribution f k is clearly equal to the equilibrium product measure 

h ~ ~ k\N-ky (3 - 10) 

When the number N of particles is large, the above distribution exhibits a segregation 
phenomenon for any positive value of b (see [7] for a more detailed analysis, including 
an asymptotic analysis of the amplitude Cn)- The most probable configurations are 
those where almost all the particles are at one site, i.e., either k <ti N or N — C iV. 
This simple example confirms that the target model, just as the ZRP, manifests a trend 
toward segregation at high density. It therefore suggests that the preferred scenario is 
that of an alternating structure. 

This picture can be corroborated and made more quantitative by means of the 
following dynamical stability analysis of the uniform situation. For the sake of generality, 
in this part of the analysis we deal with the partially asymmetric target process with 
bias p. Consider a configuration of the uniform scenario. All the sites have very high 
local densities (N m ) = p m (t). These local densities obey the exact rate equation 

^ = p((l - S Nm _ lfi )v N J + q((l - 5 Nm+1 ,o)v N J 

- p((l - S Nmfi )v Nm+1 ) - q((l - <$jv mi o)ujv m _i>- (3.H) 

In the high-density regime it is legitimate to simplify the above equation in several 
respects. The probability that a site is empty is negligible, whereas (i>jv m ) ~ v^ m ) 
1 — bj (N m ). We are thus left with 

^«6(-P- + -?--J-y (3.12) 

dt \Pm+l Pm-1 PmJ 

Let us furthermore assume that the density profile is close to being constant, i.e., 

Pm = P + 0~Pm, (3.13) 

with 5p m p. The rate equation (3.11) can then be linearized as 

d 8p m ^ b ^ _ p § pm+1 _ q 5 Pm ^) , (3.14) 
dr p z 

The component of the density profile at wavevector K therefore grows exponentially 
in time as Sp(K) ~ exp(a(K) t), where the characteristic rate a{K) is given by the 
dispersion relation 

a(K) = j 2 (l-pe- iK -qe iK ). (3.15) 
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The growth rate s(K) is given by the real part of cr(K), which reads 

s(K) = - 2 sin 2 f, (3.16) 
p A 2 

irrespective of the bias p. The expression (3.16) for the growth rate of fluctuations 
around the uniform situation is manifestly positive for all values of the wavevector K. 

The uniform scenario is therefore fully linearly unstable. As a consequence, the 
alternating scenario is the preferred one. Furthermore, the alternating structure is 
already appearing as the most favored one within the stability analysis. The most 
unstable mode indeed corresponds to K — tt, i.e., an alternating density modulation 
of the form 5p m ~ (— l) m . A uniform high-density initial configuration is therefore 
expected to smoothly relax to an alternating one by the dynamics of the asymmetric 
target process. This observation deserves, however, to be complemented with the 
following caveat. Consider the symmetric target process, corresponding to p — 1/2. 
The stability analysis still leads to the expression (3.16), so that the alternating mode 
is still the most favored one. On the other hand, the stationary state is known to 
be described by a product measure, and especially to have a single condensate. This 
suggests that the relaxation dynamics of a uniform high-density initial configuration 
will exhibit two stages: first, a rather fast relaxation to an intermediate structure with 
alternating fluctuations, then, a coarsening evolution of the system toward its true fate, 
by the merging of the excess particles into fewer and fewer condensate precursors. This 
two-stage relaxation should also hold for the dual ZRP with a uniform high-density 
initial configuration. The dynamical stability analysis of the uniform situation in the 
ZRP indeed yields exactly the same expression (3.16), again irrespective of the bias p. 

To close up, we mention that the expression (3.16) also yields some hints on the 
time scales involved in the model. The characteristic time of the most unstable mode, 
T\ oc = l/ s ( 7r ) ) gives an estimate of the local relaxation time of the alternating structure, 
at the spatial scale of two consecutive sites. On the other hand, the long-distance 
behavior of the dynamics is antidiffusive. We have indeed formally s(K) « —DK 2 as 
K — > 0, with a small negative diffusion coefficient D = —b/(2p 2 ). For a large but finite 
system made of M sites, with periodic boundary conditions, the longest characteristic 
time, T g i = 1/s(2tt/M), gives an estimate of the global relaxation time of the structure 
as a whole. In the high-density regime, the characteristic times thus defined scale as 

Tioc « 4> T gl « (3.17) 
2b g 2n 2 b v ; 

The predicted divergence of both characteristic times with density provides an a 

posteriori confirmation that the high-density regime is indeed the most interesting one. 



4. One-dimensional target process: numerical results 

In this section we complement our analysis of the target process in the asymmetric 
one-dimensional case, by means of numerical simulations and scaling arguments. 
For definiteness we restrict the study to the totally asymmetric situation {p = 1). 
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Furthermore, we set once for all 6 = 4. We successively consider features of the transient 
dynamics, of the stationary-state measure, and of the stationary-state dynamics. The 
main focus is on the scaling behavior of quantities of interest in the high-density regime. 
Two types of initial conditions are considered: 

• Deterministic initial condition: the occupations of all the sites are set equal to 
N m = P (provided the density p is an integer). 

• Random initial condition: the occupations N m are drawn independently at random 
from the geometric distribution (2.35) at density p. 

We start by investigating the early stage of the dynamics, where relatively fast 
rearrangements of particles bring the system to a locally stationary state. This stage of 
the relaxation dynamics can be monitored by means of a local probe at one site. We 
choose the reduced second moment of the occupations: 

*(0 = 4E* 2 /*(f), (4.i) 

P k>0 

where fk{t) is the time-dependent distribution of the occupations. The initial values of 
this quantity are K(0) = 1 for the deterministic initial condition, and K(0) = (2p + l)/p 
for the random one. 

Figure 2 shows a plot of K(t), measured by means of a numerical simulation, against 
time t, for p = 10, and a deterministic and a random initial condition. Each series of data 
are obtained by averaging over sufficiently many histories in order to obtain a smooth 
signal (10 4 histories of a system of 10 4 sites in this case). The data exhibit a rather 
fast rise from their initial values, and converge to a common limiting value, K pa 2.82, 
extracted from data for much longer times, and shown as a dashed line. The data for the 
deterministic initial condition (lower curve) increase as a monotonic function of time 
and present a rather sharp shoulder, whereas those for the random initial condition 
(upper curve) exhibit a non-monotonic behavior with a very flat maximum. 

The characteristic time of the rise observed in K(t) gives a measure of the local 
relaxation time. More precisely, we define the local time T\ oc by the condition that K(t) 
is near the middle of its rise, i.e., K{T\ oc ) = 2 for the deterministic initial condition, and 
^(^loc) — 2.5 for the random initial condition. Figure 3 shows a plot of the numerical 
values of the local time so defined, divided by p, against p, for both types of initial 
conditions. Here and in subsequent figures, statistical errors are comparable to the 
symbol size. The least-squares fits of the two series of data suggest a growth of the form 

T loc « V + Bp (4.2) 

in the high-density regime of interest. The amplitudes A and B depend on the initial 
condition. The least-squares fits shown on the plot yield A pa 0.095 for a deterministic 
initial condition and A pa 0.039 for a random initial condition. These numbers, and 
especially the first one, are comparable to the rough estimate coming from (3.17), i.e., 
A = 1/(26) = 1/8 = 0.125. 
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Figure 2. Plot of the reduced second moment K(t) of the occupation 
distribution in the one-dimensional fully asymmetric target process against 
time t, for 6 = 4 and p = 10. Lower curve: deterministic initial condition. 
Upper curve: random initial condition. Horizontal dashed line: common 
limiting value K m 2.82. 
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Figure 3. Plot of the local relaxation time T\ oc , divided by p, against density p. 
Upper data (empty symbols): deterministic initial condition. Lower data (full 
symbols): random initial condition. Straight lines: least-squares fits with 
respective slopes 0.095 and 0.039. 

The initial condition is irrelevant for what concerns later stages of the dynamics, 
which correspond to the emergence of global features of the stationary state. Hereafter 
we choose to work with a random initial condition. Before we turn to an analysis of 
these late stages, it is worth taking a glance at the spatial structure of the stationary 
state in the high-density regime. 

Figure 4 shows a typical stationary occupation profile for p = 50. In order to 
better reveal the alternating structure predicted in Section 3, we have plotted (— l) m N m 
against the position m of the site. The alternating structure clearly emerges from 
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Figure 4. Plot of ( — l) m N m against the position m of the site, for a 
typical stationary-state configuration for p = 50, emphasizing the alternating 
structure. 



this representation: domains where the signal is positive (resp. negative) correspond to 
domains where the B-sites are the even (resp. odd) sites. 

In order to turn this observation into a quantitative measurement, let us introduce 
the concept of defects. Roughly speaking, a defect is a site around which the structure 
is not perfectly of the alternating form BSBSBSBS... and a domain is any stretch 
between two consecutive defects. More precisely, the site m is considered as a defect 
whenever its occupation N m is neither a maximum nor a minimum of the density profile. 
Equivalently, N m is between N m _ 1 and N m+1 , i.e., the product (N m — N m _i)(N m — N m+1 ) 
is negative. This definition pinpoints 24 defects in the configuration shown in Figure 4. 
This number is slightly above the number of domains visible with the naked eye, i.e., 18, 
because some of the domains are microscopic. Most defects can be viewed either as a 
BB sequence or as an SS sequence. We shall return later on to the dynamics of these 
defects. The density of defects R, i.e., the mean number of defects per unit length, reads 

R = (G(-(N m - N m ^)(N m - N m+1 ))), (4.3) 

where the Heaviside step function on the integers is defined as 

. . f 1 for n > 0, , AA . 
e( "» = {o for»<0. (44) 

The inverse of the density of defects, 

{ = i, (4-5) 

is interpreted as the mean size of a domain or, equivalently, as the coherence length of 
the alternating structure. 

Let us now return to the late stages of the dynamics, starting from a random 
initial condition. In order to characterize the growth of the alternating structure, we 
have measured the time dependence of the mean domain size £(£) = 1/R(t). Figure 5 
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Figure 5. Plot of the mean domain size £(t) against the reduced time (t/T\ oc ) 1 / 2 , 
for various values of the density p, indicated on the curves. The dashed straight 
line starting from £(0) = 3 illustrates the coarsening law (4.6). 



shows a plot of £(£), for various values of the density p. For a random initial condition, 
we have £(0) = 3. Consider indeed the initial values of iV m _i, N m , and N m+ i. The 
probability that these three independent random numbers obey either of the inequalities 
N m -i < N m < N m+ i or N m -i > N m > N m+ \ is equal to 1/6 in the high-density regime 
(neglecting the fact that these are integer variables, which may coincide with a small but 
nonzero probability). Hence R(0) = 2 x 1/6 = 1/3 and £(0) = l/R(0) = 3. The data 
are plotted against the reduced time variable (t/T\ oc ) l l 2 , for each value of the density p, 
where T\ oc is taken from the lower data of Figure 3. The observed common initial linear 
behavior, shown as a dashed straight line starting from the known value £(0) = 3, 
demonstrates that the mean domain size grows according to the coarsening law 

m ~ {t/T loc f/\ (4.6) 

before it saturates to a density-dependent stationary-state value, simply denoted by £. 
The duration of the coarsening process, before the stationary state is reached, defines 
the global relaxation time T g \ of the problem. By inverting the relation (4.6), we predict 
that the latter time grows as T g \ ~ Ti oc £ 2 . This relation between T loc and T g \ is in 
agreement with (3.17), where the stationary-state mean domain size £ plays the role of 
the system size M. 

We now investigate a few characteristic features of the nonequilibrium stationary 
state of the asymmetric target process, emphasizing their scaling behavior at high 
density. We start with the mean domain size £. Its scaling behavior at high density can 
be predicted by the following argument. Roughly speaking, a defect can be thought 
of as a B-site whose occupation is accidentally as small as that of an S-site, i.e., 
finite. Anticipating the scaling law (4.9), and assuming a linear rise for the scaling 
function F(x), we find that the probability of such an event scales as 1/p 2 . This 
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Figure 6. 
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Plot of the mean domain size £ in the stationary state against 
Full line: second-degree polynomial fit to the data with leading 
3 = 416 x 10~ 6 . 



argument leads to an asymptotic quadratic growth of the mean domain size of the form 



Figure 6 shows a plot of the stationary-state value of £ against density p, for 
densities up to p = 200. The second-degree polynomial fit to the data is compatible 
with our expectation (4.7). The numerical value for the prefactor, S as 4 x 10~ 4 , is 
however found to be very small. Partly as a consequence of this smallness, the data 
exhibit large corrections to the above asymptotic law for values of the density accessible 
to numerical simulations. At variance with the case of 7i oc , plotted in Figure 3, we found 
no way to unambiguously characterize these corrections. The quadratic growth (4.7) 
of £ corresponds to a very fast growth of the global relaxation time: 



Figure 7 shows a plot of the stationary-state values of the current J and of the 
reduced second moment K of the occupation distribution, against the stationary-state 
density of defects R = l/£. The second-degree polynomial fit for J yields Joo ~ 0.250, 
in excellent quantitative agreement with the limiting value J = 1/6 = 1/4 predicted 
in (3.9). A third-degree polynomial fit for K yields the limiting value m 2.98. This 
number is very close to the value = 3 corresponding to the trial scaling function 
Fq(x) given in (4.11). This proximity suggests that the value = 3 could be exact. 
Imposing the constraint = 3 indeed hardly alters the quality of the fit (dashed line). 

In complete agreement with the alternating scenario depicted in Section 3, the 
occupation distribution fa in the stationary state consists of two distinct components 
with equal weights, which respectively describe S-sites and B-sites. These components 
are emphasized in Figure 8, where the data for fk for the same values of the density, 
p = 50 and 100, are plotted in two different ways. The left panel shows a plot of In fa 
against ln(A; + 1) for moderate values of the occupation k, up to 20. These are, roughly 




(4.7) 



T g \ ~ Ti oc £ 2 ~ p 6 . 



(4.8) 
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Figure 7. Plot of stationary values of the current J (left) and of the reduced 
second moment K (right), against the stationary density of defects R = l/£. 
Symbols: data for densities ranging from p = 30 to 200. Full lines: polynomial 
fits, of degree 2 for J and 3 for K, yielding the extrapolated values 0.250 
and Koo « 2.98. Dashed line: constrained polynomial fit of degree 3 for K, 
imposing = 3 (see text). 




Figure 8. Stationary-state occupation distribution The left panel shows 
data for k <C p (S-sites). Empty symbols: plot of ln/fc against In (A; + 1) for 
p = 50 (upper data) and 100 (lower data). Full symbols (labeled ZRP): plot of 
ln(/fc/2), where fk is the critical distribution of the dual ZRP. The right panel 
shows data for k comparable to the density p (B-sites). Full lines: scaling 
plots of the product pfk against k/p for p = 50 (upper data) and 100 (lower 
data). Symbols: minimum value of fk- Dashed line: trial scaling function Fq(x) 
defined in (4.11). 



speaking, the S-sites. The data for the smaller values of k are very close to 1/2 times 
the critical occupation distribution of the dual ZRP or, equivalently, of the symmetric 
target process, given by (2.41), also shown on the plot (full symbols). The range of 
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values of k over which the agreement holds is observed to get larger for larger densities. 
This is a convincing confirmation of the prediction made in Section 3 that the S-sites 
are critical in the stationary state. The right panel shows a scaling plot of the product 
pfk against k/p, for larger values of k, comparable to the density p. These are, roughly 
speaking, the B-sites. The data exhibit a scaling law of the form 

,pj ■ (49) 

The scaling function F(x) obeys the sum rules 

rOO POO POO 

/ F(x)dx = l/2, / xF(x)dx = l, / x 2 F{x) dx = K^. (4.10) 
Jo Jo Jo 

The first two equalities express that the fraction of B-sites is 1/2 and that their mean 
occupation is 2p, whereas the third one is a rewriting of the definition of K^. The 
scaling function is observed to be rather uniformly well approximated by the trial scaling 
function 

F (x) = \xe~ x , (4.11) 

shown on the right panel of Figure 8 as a dashed line. The scaling function (4.11) 
corresponds to = 3. The closeness of this number to the extrapolated value 
i^oo ~ 2.98 opens up the possibility that the scaling function F(x) is exactly given 
by F (x). In any case, the linear rise of the trial function F (x) seems to be shared by 
the true scaling function F(x). 

Both components of the occupation distribution shown in Figure 8, respectively 
corresponding to S-sites and B-sites, are separated by a minimum in the occupation 
distribution. In the high-density regime this minimum takes place for a crossover 
occupation k+ such that the estimates l/k b (see (2.43)) and k/p 2 (assuming a linear 
rise for the scaling function F(x)) are comparable. We thus obtain 

fc*~p 2/(b+1) , h,~ P - 2b/(b+1) . (4.12) 
These estimates make sense as soon as the crossover occupation obeys k± <C p. We 
thus recover the condition b > 1 for the validity of the alternating scenario. The mean 
occupation p$ of the S-sites has the finite asymptotic value p c for b > 2, whereas it 
scales as ps ~ k 2 ~ b ~ p 2 ( 2 ^ 6 )/( 6+1 ) in the high-density regime for 1 < b < 2. 

We close up this section with an investigation of the stationary-state dynamics of 
defects. For a large but finite density p, there is a small density of defects i?=l/£~l/p 2 
(see (4.7)). These defects cannot stay immobile. 

• Consider indeed a BB defect, made of two consecutive B-sites. The current between 
the two B-sites, J BB = 1, exceeds the mean current J = 1/b through the system. 
As a consequence, particles flow from the left B-site of the defect into the right one 
at a rate uj bb = J BB — J = (b — l)/b. After a time of the order of 

TBB „ 2L = »L, (4.13) 
u BB 6-1 

the left B-site is emptied. This is the first reaction of (4.14). An SS defect is thus 
formed one site to the left of the original BB defect. 
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Figure 9. Space-time plot of the stationary domain pattern for p = 100. Black 
(resp. white) areas show the regions of space-time where the S-sites are the odd 
(resp. the even) sites. The slope of the straight line in the right part of the plot 
yields V » -0.0072. 



• Consider now an SS defect, made of two consecutive S-sites. The current between 
the two S-sites, Jss ~ ^-/b 2 , is smaller than the mean current J. As a consequence, 
particles flow from the B-site to the left of the defect into the left S-site at a rate 
Uss — J — Jss ~ (b ~ l)/^ 2 - The left S-site is thus soon (i.e., after a time which 
does not grow proportionally to p) turned to a B-site. This is the second reaction 
of (4.14). A BB defect is thus formed one site to the left of the original SS defect. 

The discussion can be summarized in the form of the following reactions 

BB -> SB, SS -> BS. (4.14) 

The typical history of a single defect therefore looks as follows, where time runs from 
bottom to top, for the sake of consistency with Figure 9: 

BSBSBBSBSBSBSBSBSBSB 

BSBSBSSBSBSBSBSBSBSB 

BSBSBSBBSBSBSBSBSBSB 

BSBSBSBSSBSBSBSBSBSB 
t T BSBSBSBSBBSBSBSBSBSB 

BSBSBSBSBSSBSBSBSBSB 

BSBSBSBSBSBBSBSBSBSB 

BSBSBSBSBSBSSBSBSBSB 
n — > 

The time it takes for a defect to move two sites to the left is therefore equal to r BB 
on average. As a consequence, the pattern of defects and domains is advected with an 
upstream (negative) velocity V ~ — 2/t B b, i-e., 
b-1 



V 



bp 



(4.15) 
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This scenario is confirmed by Figure 9, showing a space-time plot of the stationary 
dynamics of the domain pattern for p = 100. Equation (4.15) predicts V ~ —0.0075, in 
reasonably good agreement with the observed value V ~ —0.0072. The plot shows that 
the whole advected domain pattern behaves more or less as a rigid body over spatial 
scales much larger than the mean domain size £. 

5. Two-dimensional target process 

We now consider the canonical target process defined by the rate (2.33), on the square 
lattice with unit vectors e 1; e 2 . In order to have a genuine two-dimensional model 
and to maximize the asymmetry, we choose the following rule: particles hop either 
East (displacement ei) or North (displacement e 2 ) with equal probabilities. In other 
words, if the departure site is d = (to, n), the arrival site is chosen to be either 
a = (to + l,n) = d + ei, or a = (to, n + 1) = d + e 2 , with probability 1/2. The 
bias, i.e., the mean displacement proposed to a particle, b = (e x + e 2 )/2, is along the 
North-East direction. 

5.1. Heuristic argument 

In the one-dimensional situation, the existence of a conserved current was instrumental 
in order to discriminate between possible scenarios for the stationary state of the model. 
In the present situation, however, the current J is a two-dimensional vector. The 
condition that J be conserved in the stationary state is less stringent than in the 
one-dimensional case. On spatial scales much larger than the lattice spacing, it is 
reasonable to use the continuum formalism. Within this framework, the conservation 
law reads V • J = 0. The symmetry of the dynamical rules implies that J is aligned 
with the bias b. Current lines are therefore parallel straight lines along this direction. 
The conservation law implies that the magnitude of the current is constant along each 
current line, but may well vary in the transversal direction from one current line to 
another. In particular the existence of an extended condensate is allowed, in the form 
of a one-dimensional structure elongated along the direction of the bias. This opens 
up the possibility of having an unconventional type of condensation transition. Such a 
phenomenon is indeed observed in the numerical simulations described hereafter. 

5.2. Numerical results 

We choose once for all the value b = 4 in the numerical simulations of the two- 
dimensional target process defined above. 

In analogy with the one-dimensional case, we begin with the transient dynamics, 
starting from a random initial condition. We monitor the local relaxation by means of 
the reduced second moment K(t) of the occupation distribution, introduced in (4.1). 
Figure 10 shows a plot of K{t) against time t. For all values of the density p, K(t) is an 
increasing function of time, starting from its initial value K(0) = (2p+ I)/ p. The most 
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Figure 10. Plot of the reduced second moment K(t) of the occupation 
distribution in the two-dimensional target process against time t, for various 
values of the density p, indicated on the curves. Left: data for moderate 
values of time. Right: data for longer values of time, multiplied by density p. 
The parallel straight dashed lines, meant as a guide to the eye, have a slope 
9.1 x 1(T 3 . 



significant feature to be observed on the data is the following: K(t) saturates to a finite 
limiting value, denoted K, for values of the density up to the threshold value p 6. 
We have K ps 4.8 for p = p . To the contrary, K(t) grows indefinitely for larger values 
of the density. The right panel of Figure 10 demonstrates that the asymptotic growth 
law of K(t) for all p > p is of the form 

(It 

Kit) » — , (5.1) 
P 

where C^9x 1CT 3 is taken from the slope of the parallel dashed lines. 

Pursuing along the lines of our investigation of the one-dimensional case, we 
introduce a local relaxation time T\ oc , denned by the condition K(T\ oc ) = 4. Figure 11 
shows a plot of T\ oc so defined, divided by p, against p. The least-squares fit suggests a 
growth of the form (4.2), with A m 0.63, in a wide range of values of the density. Notice 
that the local relaxation time exhibits no visible singularity at the threshold density po- 

We now turn to a more accurate determination of the threshold density po, using 
the pre-asymptotic growth of K(t) in the intermediate time regime, where K(t) has 
already departed in a significant way from its plateau value K , but not yet reached 
the asymptotic linear growth (5.1). The right panel of Figure 10 demonstrates that 
this regime lasts longer and longer as the threshold density po is approached. This 
observation is turned to a quantitative measurement by defining the global relaxation 
time T g i by the condition K(T g i) = Kq + AK, where we set Kq = 4.8, whereas the 
choice AK = 50/ p incorporates the form (5.1) of the asymptotic growth law. Figure 12 
shows a plot of the reciprocal of the global time T g \ thus defined, against density. 
The data convincingly demonstrate that the global time diverges at some non-trivial 
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Figure 11. Plot of the local relaxation time Tj oc , divided by p, against density p. 
Full straight line: least-squares fit with slope A « 0.63. 




Figure 12. Plot of the reciprocal of the global time T g \ against density p. Full 
line: second-degree polynomial fit yielding po = 6.0 ± 0.1. 



threshold density po- A crossover of the data toward another type of asymptotic behavior 
indeed seems extremely improbable, in view of the accuracy of the available data. The 
second-degree polynomial fit to the data shown on the plot provides a rather accurate 
determination of the threshold density, 

p = 6.0±0.1, (5.2) 

as well as an evidence that the global time diverges linearly as the threshold density is 
approached from above, as 

Tgl w To (5.3) 
P - Po 

with a rather large prefactor T ~ 28 000. 

The threshold density po is the maximal density for which a homogeneous fluid 
phase is stable. At variance with the critical density p c of the symmetric target process 




Figure 13. Logarithmic plot of occupation distributions. Upper curve: fully 
asymmetric target process in two dimensions at its threshold density pQ. The 
dashed line, meant as a guide to the eye, has a slope — p w —0.06. Lower curve: 
critical occupation distribution of the dual ZRP, given by (2.41). 



and of the dual ZRP, p rather appears as a dynamical threshold. This viewpoint is 
corroborated by the fact that the stationary state at density po does not exhibit any 
critical feature. The distribution of the site occupations at the threshold density, shown 
in Figure 13, has an exponential fall-off of the form ~ exp(— pk), with p pa 0.06, at 
least in the accessible range of values of the occupation. We checked that the data are 
not affected in an appreciable way by finite-size effects in the range considered. The 
critical occupation distribution (2.41) of the dual ZRP is shown on the same plot as a 
comparison. The latter distribution has a much smaller density p c = 1/2, some 12 times 
smaller than the observed threshold density (5.2) of the asymmetric two-dimensional 
model, but a slower power-law fall-off, so that the distributions eventually cross each 
other. 

The progressive emergence of highly occupied coherent structures, which are 
strongly elongated along the direction of the bias, is illustrated in Figure 14, showing 
two snapshots of the coarsening regime of a sample of size 100 x 100 at density p = 20. 
The filled symbols show the 1% most occupied sites. The visible structures clearly are 
precursors of the extended condensate mentioned above. The background density of 
the fluid phase besides these structures is found to be much smaller than the threshold 
density p , and comparable to the critical density p c — 1/2 of the dual ZRP. We shall 
return to this point in more detail below (see Figure 16). 

We now study the typical characteristic sizes (width and height) of the extended 
condensate in the stationary state of a finite system, and of its precursors in the 
coarsening regime of an infinite system. The width W of a condensate is denned as 
the number of sites which take part in the condensate, whereas its height H is the mean 
number of particles per site in the condensate, so that the product HW gives a measure 
of the number of particles involved in the condensate. We first consider the coarsening 
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Figure 14. Plots of the 100 (i.e., 1%) most occupied sites of a sample of size 
100 x 100 at density p = 20 in the coarsening regime. Left: t = 5 x 10 4 . Right: 
t = 10 5 . 



regime of an infinite system. The condensate precursors shown in Figure 14 are expected 
to be characterized by a typical width W(t) and height H(t), with both scales growing 
with time. The contribution of these precursors to the reduced second moment of the 
occupation distribution can be checked to scale as K(t) ~ H(t)/p, irrespective of the 
width W(t). The growth law (5.1) therefore implies that the height of condensate 
precursors grows linearly in time, according to 

H{t) w Ct, (5.4) 

with C ~ 9 x 10 -3 . This asymptotic coarsening law is expected to hold for any 
density p > p$. The behavior of the width W of condensate precursors can only 
be investigated in an indirect way, by means of finite-size scaling. We therefore 
consider finite systems, namely square samples of linear size L, with periodic boundary 
conditions. Figure 15 shows numerical data concerning the stationary state of finite 
systems against their linear size L, at fixed density p = 20, well above the threshold 
density p . The left panel shows the stationary-state value K L of the reduced second 
moment of the occupation distribution. The right panel shows the characteristic 
relaxation time T L) defined by the condition K(T L ) = (K L +K )/2, again with K = 4.8. 
The data for both quantities clearly exhibit a linear growth with the size L. The first 
of these growth laws implies H L pa pK L ~ L. This is in accord with the expectation 
that typical stationary-state configurations have a single and roughly system-spanning 
extended condensate, for which Hl ~ L. Furthermore, as the number of particles 
involved in the condensate scales as HlWl ~ L 2 , we have Hl ~ Wl ~ L. Let us 
now make the finite-size scaling assumption that H(t) and W(t) become respectively 
comparable to H L and W L for a time t comparable to the relaxation time T L . This 
yields the scaling law H(t) ~ t, already known (see (5.4)), and the prediction W(t) ~ t. 
Finally, we have also measured the background density p L g of the fluid phase in the 
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Figure 15. Plots of data concerning the stationary state of finite systems at 
fixed density p = 20, against their linear size L. Left: stationary-state value Kl 
of the reduced second moment of the occupation distribution. Right: relaxation 
time Tl- Full straight lines: least-squares fits with respective slopes 1.79 and 
2 660. 



stationary state of finite samples at density p = 20. This quantity is algorithmically 
denned as follows. For any intercept k — 1, . . . , L, consider the total number of particles 
in the diagonal array with intercept k, i.e., with equation n = m + k (mod. L): 

L 

Nk=J2 N m,rn+k- (5.5) 
m=l 

The largest of these L numbers, A/" max , corresponds to the diagonal array occupied by 
the extended condensate. It is overwhelmingly larger than the others, as it scales as 
A/" max ~ L 2 . The other (L — 1) numbers Af k represent the fluid phase, and therefore scale 
as PiL. We are thus naturally led to define the background density as 

bg _ N - Mnax ( r fi x 

Pl ~ L(L - 1) ' (5 - 6) 
where N is the total number of particles in the system. Figure 16 shows a plot of the 
stationary-state background density p^ g for finite samples of size L, against 1/L. The 
data demonstrate that is smaller than unity, and therefore much smaller than the 
mean density p = 20, as soon as L > 7. The segregation phenomenon is therefore 
already fully at work for rather small system sizes. From a quantitative viewpoint, 
a second-degree polynomial fit to the data yields the extrapolated value p bg ~ 0.48. 
The limiting value thus obtained is remarkable close to the critical density p c — 1/2 of 
the dual ZRP. Imposing the constraint p bg = p c = 1/2 indeed hardly changes the fit 
(dashed line). This agreement strongly suggests that the background fluid phase of the 
two-dimensional target process above its threshold density is characterized by the critical 
occupation distribution of the dual ZRP, just as the S-sites of the one-dimensional case 
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Figure 16. Plot of the background density p L g in the stationary state of finite 
systems at fixed density p = 20, against their reciprocal linear size 1/L. Full 
line: second-degree polynomial fit yielding the extrapolated value p° g « 0.48. 
Dashed line: constrained polynomial fit of degree 3 imposing p hg = p c = 1/2. 



in the high-density limit. As a consequence, the density p hg of the fluid phase: 



[p c for p > po, 
has a discontinuous jump at the threshold density p . 

6. Discussion 

In this work we introduced a novel example of a migration process, the target process. 
We then studied in detail the structure of the nonequilibrium stationary state of 
the asymmetric target process, the main focus being on the fate of the condensation 
phenomenon. 

The stationary-state measures of migration processes do not have a product form in 
general. The symmetric target process, though, has the same stationary-state product 
measure as the corresponding dual ZRP. In particular, the so-called canonical target 
process, defined by the rate (2.33) dual to the ZRP with rate (2.32), has a continuous 
condensation transition at a finite critical density p c whenever b > 2, with a macroscopic 
condensate occupying a single site for densities p > p c , irrespective of the dimensionality 
of the system. The asymmetric target process has a fluctuating stationary state with 
non-trivial spatial and temporal correlations, whose qualitative features depend on the 
dimensionality. Our main effort in the present work consisted in characterizing this 
nonequilibrium stationary state, including its dependence on the dimensionality. 

We have reached a complete understanding of the one-dimensional target process, 
especially in the high-density regime of most interest. We showed, by exploiting the 
existence of a conserved current, that the asymmetric canonical target process has no 




for p < po, 



(5.7) 
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condensation transition, and remains homogeneous at any finite density. In the high- 
density regime, an alternating scenario prevails for b > 1: typical configurations consist 
of long alternating sequences BSBSBSBS... of highly occupied B-sites, and less occupied 
S-sites, whose occupation distribution coincides with the critical distribution of the dual 
ZRP. The coherence length (mean domain size) of this alternating structure diverges as 
£ ~ p 2 . We also gave a characterization of the scaling behavior of many other quantities 
in the vicinity of the 'infinite-density fixed point'. 

For the asymmetric target process in higher dimensions, we argued that the 
condensate must be extended and have the form of a one- dimensional structure elongated 
along the direction of the bias. In the two-dimensional case, numerical simulations 
performed for 6 = 4 show that the model exhibits an unconventional condensation 
transition at the density po ~ 6. This density, which is much larger than the critical 
density p c — 1/2 of the dual ZRP, appears as a dynamical threshold: it is the maximal 
density at which a homogeneous fluid phase is dynamically stable. For p > p , the 
predicted extended condensate is observed, whereas the background fluid phase again 
appears as critical. This picture seems to be generic for higher-dimensional systems. 
Preliminary numerical simulations of the asymmetric target process on the three- 
dimensional cubic lattice (where the displacement is along either of the unit vectors 
ei, e 2 , e 3 with equal probabilities), indeed show that the overall picture is quite similar 
to the two-dimensional one. The global relaxation time T g i is again found to diverge 
according to (5.3), with p ~ 34 for 6 = 4. 

Let us mention that another mechanism leading to an extended condensate in 
a class of mass transport models in one dimension has been reported recently [20]. 
There, the nonequilibrium stationary-state measure is a product whose factors involve 
the occupations of two consecutive sites. In the condensed phase, those models exhibit 
an extended condensate, whose height and width scale as H N ~ W N ~ N 1 / 2 for a finite 
system of N sites. These scaling laws are formally identical to those found in the present 
work. 

The existence of a threshold density po at which the background density has a 
discontinuous jump (see (5.7)) is reminiscent of what occurs in the model studied in [21], 
namely a ZRP with two species of particles, and with rates such that the stationary- 
state measure does not have a product form. When the densities p^ and p^ of the 
two species are equal, the behavior of the system is qualitatively the same as that of 
the canonical ZRP (with one species). In particular the system has a continuous phase 
transition at some critical density p c , from a fluid phase to a condensed phase with 
critical background. The general situation where the two densities are different however 
drastically departs from this known scenario. If either of the two densities (p^\ say) 
is kept fixed at a value larger than p c , on increasing the other density p^ 2 \ the system 
remains homogeneous as long as p^ is less than a threshold value p^ which depends 
on p^\ At this threshold the system undergoes a discontinuous transition from an 
imbalanced fluid phase, where both species have densities p^ and p^ 1 larger than the 
critical density, to an imbalanced condensate coexisting with a balanced critical fluid 
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with densities p^ = p^ = p c . 

Finally, the observed rapid growth of the threshold density po with the dimen- 
sionality raises the question of the behavior of the asymmetric target process in high 
dimensions. At this point let us emphasize that the absence of a stationary-state product 
measure for the target process is a rather subtle effect which needs the conjunction of 
several ingredients, and chiefly the presence of a bias. This feature cannot be present 
in mean-field geometries such as the complete graph, so that the dynamical threshold 
behavior of the model in high dimensions is not expected to smoothly converge to a 
well-defined mean-field limit. 
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